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Abstract 

Fermion-antifermion pair-production in the presence of classical fields is described based on the 
retarded and advanced fermion propagators. They are obtained by solving the equation of motion 
for the Dirac Green's functions with the respective boundary conditions to all orders in the field. 
Subsequently, various approximation schemes fit for different field configurations are explained. 
This includes longitudinally boost-invariant forms. Those occur frequently in the description of 
ultrarelativistic heavy-ion collisions in the semiclassical limit. As a next step, the gauge invariance 
of the expression for the expectation value of the number of produced fermion-antifermion pairs as 
a functional of said propagators is investigated in detail. Finally, the calculations are carried out 
for a longitudinally boost-invariant model-field, taking care of the last issue, especially. 

PACS numbers: 11.15.Kc, 12.38.Mh, 25.75.Dw 
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I. INTRODUCTION 



In |l| the inclusive spectra of fermion-antifermion pairs produced by vacuum polarisation 
in the presence of time-dependent classical fields were studied based on the exact retarded 
propagator. Here the generalisation to arbitrarily space-time dependent fields is presented. 
Processes of this kind in quantum electrodynamics (QED) 0, 01 as well as in quantum 
chromodynamics (QCD) are of importance for the physics of the early universe and 
strong laser fields p as well as of ultrarelativistic heavy-ion collisions and the quark-gluon 
plasma (QGP) 0. The strong classical fields are the common feature of these phenomena, 
while the detailed characteristics of the field configurations differ. Although this paper 
ultimately concentrates on ultrarelativistic heavy-ion collisions, various other kinematical 
regimes are treated as well (see section [TTJ) . 

A lot of effort is made to study the QGP's production and equilibration (if in nuclear 
collision experiments at the Relativistic Heavy-Ion Collider (RHIC) and the Large Hadron 
Collider (LHC), currently under construction. It is a widely used hypothesis that the initial 
state in heavy-ion collisions is dominated by gluons which on account of their large occu- 
pation number can be treated as a classical background field. The larger the occupation 
number Nk of the bosonic sector at a given momentum k of a physical system, i.e., the better 
Nk 3> 1 is satisfied, the more accurately it can be described by a classical field. In a heavy- 
ion collision at RHIC with = 130GeV the initial occupation number for gluons at their 
average transverse momentum |/c^| « lGeV in the center of the collision is approximately 
equal to 1.5 9]. Even if this number is not much larger than unity, the classical field as the 
expectation value of the gauge field still represents the parametrically leading contribution. 
At this scale, keeping only classical bosons is a first approximation but quantum fluctuations 
have to be investigated subsequently. 

One concrete model based on this idea is due to McLerran and Venugopalan In 
that approach, classically interacting colour-charge distributions on the two branches of the 
light-cone represent the two colliding nuclei. However, there, it turns out that quantum 
corrections become important due to kinematical factors. Nevertheless, this approach can 
be extended farther by including quantum evolution in the distribution of the classically 
interacting charges. This extension leads to the so-called JIMWLK equation [Io| . 

Be that as it may, the magnitude of the high occupation-number bosonic fields is such 
that multiple couplings to the classical field, i.e., higher powers of the product of the cou- 
pling constant g and the gauge field strength A, are not parametrically suppressed. In the 
absence of other scales they have to be addressed to all orders. Under the prerequisite of 
weak coupling the leading quantum processes concern terms in the classical action of sec- 
ond order in the quantum fields (fermions, antifermions, and bosonic quantum-fluctuations). 
These terms' coefficients give the inverse of the respective two-point Green's functions. Their 
inversion for selected boundary conditions yields the propagators of the corresponding par- 
ticles to all orders in the classical field. As an alternative, the full propagator can also be 
derived by resumming all terms of the perturbative series or by adding up a complete set of 
wave-function solutions of the equations of motion. 

Since exact solutions are not easily obtainable some general approximation schemes have 
been devised. Expansion of the propagators in powers of the field leads to the perturbative 
series which is based on the free Green's functions. Neglecting the spatial part of the field 
and the derivatives yields what is called the static approximation. Other approaches for the 
fermionic propagator can be found in [rj| . 
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In this paper the phenomenon of particle production is to be investigated based on the 
propagator in the classical field. The required correlator is obtained by direct solution 
of the equation of motion. In QED as well as QCD the direct production of fermion- 
antifermion pairs occurs (electron-positron and quark-antiquark, respectively). Caused by 
the non-linearity of the gauge-field tensor in QCD also pairs of gluonic quantum fluctuations 
are created. There, in strong fields the two types of reactions are equally parametrically 
favoured. 

In a given situation the bosonic sector could be covered sufficiently well by the concept 
of a classical field and the hard sector be accessible perturbatively |l2|. However, due to 
Pauli's principle, it exists no such concept for the fermions. Perturbative investigations can 
only describe the hard not the soft part. Hence, this paper deals predominantly with the 
description of fermions and antifermions. 

In constant fields, the Schwinger mechanism gives the exact answer to the problem of 
particle production. Frequently, its application is extended to slowly varying fields and/or 
low-energy particles. It, like most of the other schemes, is based on assumptions on time 
and/or energy scales of the situation under investigation. However, if such scales are to be 
obtained by means of a calculation it is likely that they are strongly biased by the choice of 
the approximation . 



Finally, it is necessary to decide whether expectation values are to be obtained in the 
framework of the m-m-formalism of quantum field theory - i.e., based on the advanced and 
retarded propagator and their relatives - or probabilities within that of the in-out- formalism 
-i.e., starting out from the Feynman propagator - [3,0. As the average number of produced 
pairs is the more natural observable as compared to the probability to produce exactly a 
given number n of pairs in a particular event, the first possibility will be pursued here. 

Therefore, section ITTlpresents the solution of the equation of motion for the retarded Dirac 
propagator to all orders in the field. Respective subsections include details on a generalised 
translation operator that occurs during the derivation as well as on the propagator in a large 
variety of expansion schemes and classes of field configurations. This includes longitudinally 
boost-invariant field configurations which can be of importance in the description of the 
initial phase of an ultrarelativistic heavy-ion collisions. Section III contains the application 
of the results of section II to the problem of particle production in general. This includes a 
thorough discussion of the gauge invariance of the result. In a subsection the formulae are 
evaluated for a longitudinally boost-invariant model-field. In the last chapter the contents 
of the paper are summarised. 

Throughout the paper the metric tensor is given by: g^ v = diag(l, —1, —1, —1), angular 
momenta are measured in units of h, and velocities in fractions of the speed of light c. From 
hereon, the coupling constant is included in the definition of the classical field: gA^ ld = A^ ew . 
The convention for Fourier transformations of one-point functions is: 

that for two-point functions: 
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II. THE RETARDED FERMION-PROPAGATOR 



If the fermionic propagators for the m-m-formalism are known to all orders in the classical 
field the generating functional of the Green's functions without radiative corrections can be 
written as: 

Z[fj,r)} = exp jz j d 4 xd 4 y fj <ri (x)G <nffa {x,y)T] ffa {y)^ , (3) 

with implicit sums over £ {+? — }• These indices refer to the sign of the imaginary 
part of the time component of x or y respectively, taken along the contour in the complex 
time-plane shown in figure Hence, the components of the matrix propagator G aia2 

correspond to: 

G ++ (x,y) = +(0 in \T${y)i>{x)\0 in ) 

= + (0 in | T-^(y)i;(x) \0 in ) 
G+-(x,y) = +(0in\i>(y)i/)(x)\0 in ) 

G_+(x, y) = - (0 m | if>{x)$(y) |0 m ) , (4) 

where T -1 stands for reverse time-ordering. The components of the matrix propagator G aia2 
can be linked to the advanced (A), retarded (R), and on-shell (S) propagator by: 



2G ai(72 = G s + <J\Ga + ct 2 Gr. (5) 

What is called the on-shell propagator is a solution of the homogeneous Dirac equation. The 
advanced and the retarded propagator solve the equation of motion of the Dirac Green's 
functions: 



fry ■ d{x) + 7 • A(x) — m] G(x, y) = 5^\x — y). (6) 

Let us consider homogeneous solutions Gn(x,y) of Dirac's equation, i.e., with the 5- 
distribution replaced by zero using a product as an ansatz: 



G H (x,y) = U(x - y )G H (x,y). (7) 
The differential equation obeyed by the functional U(xq — y$) is to be the following: 

[i"/ d (x)U(x - y )]G H {x,y) + i^dj{x)U{xQ - y )G H {x,y) = 0, (8) 

with a summation over j G {1,2,3} and where the solution for the boundary condition 
£7(0) = 1 is given by: 

U{x - y Q ) = exp[-(x - yohV^)]- (9) 
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This object corresponds to a generalised translation operator taking care of the relativity 
of events, as shall be explained in section III Al With equation (JHJ) and after subsequent 
multiplication by —iU(yo — £0)7° from the left, the remaining differential equation for the 
second factor Gh(x,v) of the homogeneous solution reads: 



d (x)G H (x,y) - iU(y - £o)7°[7 • M x ) ~ m}U(x - y )G H (x,y) = 0, 



(10) 



where the inverse of the generalised translation operator Q is given by inverting the sign of 
its argument: U~ 1 (xo — yo) = U(yo — x$). For the boundary condition Gh(xo = yo, x, y) = 1 
the last differential equation is solved by a path-ordered exponential [l|: 



G H {x,y) = Pexp 



d£oU(y - £o)7° [7 • A fa, x) - m] Ufa - y ) 



,'/o 



(11) 



Other boundary conditions would only lead to additional factors to the right of the above 
solution which would not result in independent expressions, because a homogeneous dif- 
ferential equation is investigated here. The homogeneous solution Gn{x,y) can be re- 
assembled from equations Q and (jlljl according to equation (J7j). The boundary con- 
dition for vanishing time difference xq — yo for all Dirac Green's functions is given by 
[linijEQ-^Q+o — lim x . ^y _o]G(x, y) = —ij°S^(x — y). Additionally, the retarded propagator 
Gr(x, y) is required to vanish for negative time differences: Gr(x, y) = if xo < yo; the ad- 
vanced for positive: Ga(x, y) = if Xo > yo- By virtue of the boundary condition chosen for 
equation (fTTj) the retarded propagator can be expressed as: 



iG R (x,y)-f° = +G H (x,y)5 (3 \x-y)9(xo 



the advanced as: 



(12) 



iG A (x,y)j° = -G H (x,y)S (3 \x - y)9(y - x ). 
With the help of the general resummation formula derived in pj: 

Pexp{£ ^o[5(eo) + C(eo)]} = 
= 7>expj£°d£o£(£o)} x 

r%o ( ryo *) { rio 

xPexp / d£oVexp\ dz B(z )\ Cfa)V exp / dz B(z Q/ 

wo ( r x ° 1 f f£° 1 

= Vexp / d^oVexpl dz B(zo) \Cfa)Vexp I dz B(z )\ 

Jyo U£o > \Jxq J 

xPexp / d^Bfa 



(13) 



(14) 



which is based on the group property valid for path-ordered exponentials: 
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Pexp |£° dtoBfa)} = Pexp |£° df £(fo)} * V ex P {£ ^(£0)} , (15) 

the generalised translation operators can be included into the principal path-ordered expo- 
nential: 

tG R {x,yh° = J ^ e +fk ^G H (x,y,k)6(x -y ) = 

(l ' k . e +iHs-ff) v exp | • T° [ 4 y 9 (a.) + y*. + 7 . A (£ 0; £) _ m l \ 0( Xo _ ^ 



(27T) 
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(16) 

Here, the ^-distribution has been replaced by its Fourier representation. The path-ordered 
exponential in equation (|16j) corresponds to a special mixed representation of the homo- 
geneous solution Gh- Only the spatial relative coordinates have been Fourier transformed 
to the momentum space. The absolute coordinate is kept in the coordinate space. This 
becomes clearer by noticing that, because of the 5-distribution introduced by the bound- 
ary conditions, x in the path-ordered exponential can always be replaced by y and also by 
(x + y)/2. Due to its mixed boundary conditions for positive and negative energy solutions, 
the Feynman propagator cannot be expressed in a similar way. Henceforth, unless otherwise 
stated, G(x, y) also without the index R shall stand for the retarded propagator. 

A. Generalised translation operator 

In one spatial dimension, the object in equation (jHJ) is linked to common translation 
operators exp{— td 3 } through: 

U 3 (+t) = exp{-t7 7 3 5 3 } = p + e~ m + p~ e +td \ (17) 

where 2p ± = 1 ± 7°7 3 are two projectors, whence (p^) 2 = p ± and p ± p T = 0. In three di- 
mensions, the operator (J§J) transforms a 5-distribution localised at x = y into: 

exp{-t7V^0r)}^ 3 >(£ -y) = j-[d t + l°l j d 3 (x)] ^ - ® - | ~ |^ ~ ^ + t} . (18) 

The effect on any other function can be investigated by a convolution of it with the previous 
equation. The corresponding integral has only support in the distance t from the point x 
along the light-cone, i.e., on a sphere of radius t. After making use of the one-dimensional 
^-distributions there remain two angular integrations <i$ which represent an average of the 
equidistant shifts over all directions (\z \ = 1): 

exp^V'c^)}/^) = \[dt + 7V'^(x)] J ^t[f(x + h) + f{x - k)}, (19) 

where x is the center of the coordinate system. The differential operator on the right-hand 
side of the last equation serves to assign the correct element of the Clifford algebra to every 
contribution of the integral. 
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B. Expansion schemes 



In this section, various expansion schemes for the retarded propagator are presented. 
This shall serve to establish the connection to the results in the purely time-dependent case 
PJ. Like there they are also to help to expose the characteristics of the full expression. 

1. Weak- field expansion 

In order to obtain the weak-field expansion, the general resummation formula (JUjl 
can be used to transform the path-ordered exponential in equation (fTBj) with the choice 

C(£o) = ^7°7 • A(£o). After expanding the result in powers of the gauge field, Gh = J2n Gh\ 
and Fourier transformation of the vector potential into a mixed representation one finds up 
to the second order: 

Gf(x, y, k) = exp{n°b J ^ - m](z„ - y„)} (20) 

x exp{i7°[7 J (kj + Kj) - m)(x - ^o)}«7°7 • A(£ , «) x 
xexp{z 7 [7^%-m](eo-2/o)}, (21) 

Gg>(*, y, k) = i r dfoi /* drj I >< 

Jy Jy J (2-K) 6 {l-ny 

x exp{i-f [-f j (kj + ^ + Aj) - m](x - £o)}«7°7 • Afa, k) x 
x exp{i7°[y (fej- + Xj) - m](£ - ^o)}«7°7 • ^(^o, A) x 
xexp{z7 [7 J '% -m](% -yo)}- (22) 

Subsequent orders are obtained analogously. Replacing the exponential functions in the 
higher orders by the lowest order expression in equation ([20)1. reproduces the usual form of 
the weak-field expansion. 

2. Strong-field expansion 

The expansion in the last section has been adapted for fields with A <C uj. In the opposite 
case A^> uj, for a strong field, the redistribution of the exponent of the retarded propagator 
according to equation (JT3j) must be such that C(£o) = il°[l^kj — m). Hence the retarded 
propagator becomes: 

-, r f x o ( r x o 

G H (x,y,k) = Vexp 



i / d£ Vexp < i / dz j [i"/ j dj(x) + 7 • A(z , x)\ \ x 
Jyo I J£o 

x 7°[7 : ' kj — m]Vexp ji J dzoj [i r y j dj(x) + 7 • A(z , x)\ 
x Pexpjz £ ^o7°[^^(x)+7-A(eo,x)]}, (23) 
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which subsequently can be expanded in powers of the term involving the mass and the 
three-momentum. The secondary path-ordered exponentials can again be decomposed into 
field insertions and generalised translation operators: 



Pexp U r d^y^djix) + 7 • A(£ , x)} 



Pexp 



' Vo 
rxo 



H dt U(x - ^ )7°7 ■ A(£o, x)U(£ - x ) 
Jyo 



U {x - y ) 



(24) 



Based on what was said about the general translation operator in section III A| the last 
expression shows how the behaviour of the retarded propagator depends on the values of 
the field at the relativistically retarded positions (with a light-like displacement from x). 
That equation (|2*4"j) corresponds really to an expansion in the on-shell energy u can be seen 
from the fact that: {^[j^kj — m]} 2 = u 2 . As the on-shell energy uj equals the (asymptotic) 
relativistic mass, the strong- field expansion represents one in the "inertia" as opposed to 
one in the "accelerations" which is the weak-field expansion. 



3. Gradient expansion 

For gauge fields whose temporal gradients are much larger than their spatial ones, a 
gradient expansion is obtained with the help of the general resummation formula (fT4*|) by 
choosing C = —^^df 



G H (x,y,k) = Vexp 



I d^ Vexplif dz Q ^°[y j kj - m + 7 • A(z Q , x)] j x 
Jyo I Ho > 

x z7°7 J 9j(x)'Pexp j z J dzo^l^kj — m + 7 ■ A(zq, x)] 
xPexp ji J X ° d&^Wkj - m + 7 • A(£ , x)]| . 



(25) 



A result analogous to the most general in |l| is obtained as the lowest order term of an 
expansion of the above equation in the spatial gradient. Higher orders contain powers of 
the spatial gradient but no generalised translation operators. This is consistent with the 
assumption of a slowly varying field: At the retarded positions it can be evaluated by means 
of a Taylor expansion. 



4- Combined strong-field and gradient expansion 

In the case, where the field is strong compared to the scale set by the on-shell energy 
u and above that has only relatively small spatial gradients, an additional expansion of 
equation in these yields: 
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V exp U r d^y^d 3 {x) + 7 • A(£ , x)} 



"Pexp 



yo 



.'7(1 



i V exp ^ z / (iz 7°7 • ^(^o, ^) \ «7°7 J dj(x)V exp ^ i / cfe 7°7 • ^.(^Oj 



X 



Pexpj* f X ° d^-A^x) 



(26) 



Finite orders of an expansion of this expression in the spatial gradients inserted into equation 
(|2~3|) lead to the generalisation of the strong-field approximation presented previously in [l[ . 

5. Abelian expansion 

This is the only expansion scheme not based on the resummation formula (|14jl. The 
iV-th order Abelian approximation can be obtained according to the following prescription: 
Starting out with equation (JTHJ), the group property (fTHjl is used in order to obtain a product 
of N + 1 path-ordered exponentials. There is a certain arbitrariness in choosing the inter- 
mediate points, but they must be ordered in time (see below). Subsequently, generalised 
translation operators are generated in every single factor by resumming the derivatives. 
Finally, the path ordering is neglected: 

G H (x,y,k) = 

= V J] exp L t dz U(& - 2 )7 Wkj + 7 ■ A(z , x) - m] U(z - &,)} U{£ v - f„ +1 ), 

(27) 



with xo = £o > Ci > ••• > £n > 6v+i = Vo because for the retarded propagator x > yo- V 
denotes that the factors are ordered with respect to the index v with the lowest index furthest 
to the left. 

For infinitely small interval widths (infinitely many factors) the exact result is reobtained, 
because in that case the exponential functions can be replaced by linear factors and hence 
also by path-ordered exponentials. This again gives the exact retarded propagator. However, 
the present Abelian approximation scheme converges faster than the expression with linear 
factors 0| . The above result can be interpreted as the propagation of the fermions over 
the interval with their average canonical momentum. The average comes ever 

closer to the actual value if the intervals become smaller. To the contrary, in the weak-field 
expansion, the particles are always propagated with their asymptotic kinematic momentum 
and in higher orders they only interact more and more often with the field. 



C. Different classes of field configurations 

In order to obtain further insight into the structure of the propagator, let us investigate 
the case where the gauge field depends mostly on the coordinates in the zero (time) and 
three (longitudinal) directions. Then the propagator is treated most conveniently in a mixed 
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representation as a function of said coordinates and the transverse momentum. With the 
choice B = — 7°7 3 c?3 the resummation formula (jl4j) leads to: 

/fttC — 
3 e+ifc ,( X3 - w ) 0(a . o _ yo) G( XOi y0; X3} fc3j fcr )[/3(x - yo) (28) 

with 

G(a;o, 2/0, £3, h, k T ) = "Pexp |z J d£ U 3 (x - £oh°h J kj - m + 7 • A(£ , a^t^o - ^o)| 

(29) 

and implicit summations over the transverse coordinates J G {1,2}. 

1. Gauge field with longitudinal and temporal components and gradients 

At first only the zero and three component of the gauge field are to be non-vanishing. 
Additionally, the particles are not to have a transverse mass. Hence, in the resulting equation 
also the kernel (C) involving the fields and momenta can be decomposed with the help of the 
projectors p . Contributions to the different subspaces do not mix. Therefore, all common 
translation operators meet their inverse after acting on the components of the gauge field 
and one gets: 

At ~° { M(so-^o) J 

+ p-e- lk3{xo ~ yo) Vexp\i f° d£ + A-(x + + Z+,x-)\, (30) 

[ JV2(y -x ) J 

with the definition \/2v± = Vq±v 3 and the change of variables to A = A(x + , x_). The path- 
ordered exponentials on the right-hand side are Wilson lines in which the A-t-component 
is integrated along the ^-coordinate, respectively. As a consequence of relativity, the end- 
points of the path in the ^-direction depend on the value of the x^-coordinate, respectively. 
The two contributions which come exactly from on the backward light-cone are schematically 
depicted in figure El 

Loosening the requirements on the system by admitting a non-zero transverse mass leads 
to four contributions in the exponent of the analogous expression: 



Gx T =o( x o^yo,X3,k 3 , k T ) = Pexp 



,0 

i d£ x 

"'yo-zo 

x (p+ {v / 2A+(x+, x_ + v^fo) + h} p + + p- [V2A_(x + + \/2£ , x_) - A; 3 } p~+ 
+ p + i°{i J kj - m}e +2 ^ x) p- + p- 7 °{7 J fcj - m}e~ 2 ^ x) p + )] (31) 

where use has been made of the idempotency of the projectors and of the relation 
P ± 1°{l J kj — m} = r )°{ r j J kj — m}p T . In the terms involving the transverse mass - note 



10 



that [7°(7 J /cj — m)] 2 = m^ 2 - the translation operators do not meet their inverse and hence 
do not drop out. A resummation of this expression for small transverse mass m-r yields: 



G a t =6 ( x o, Vo, x 3 , fa, k T ) = Pexp 



: / d&G^Zxixo, x + £o, x 3 , fc 3 ) x 

Jy -x a t~U 

x (p + l°{l J kj - m}e +2 Z° d3(x) p~ + p~7 { 7 J fcj - m}e- 2fo93(a V) x 

x G Z=^ X ° + £o, x , x 3) fe 3 )] G^: s (x , y , x 3 , h), (32) 

where G^ T ~5(x , y , £3, ^3) has been taken from equation (|50J). Therefore, the lowest (zero) 
order term in an expansion of the above expression in the transverse mass equals the result 
for vanishing transverse mass. The first order expression reads: 

G^(x ,y ,x 3 ,k 3 ,k T ) = £ ie =F*3(*o-i») f° d^ e^ +id ^°Gf% (33) 

where the sum over the contribution with the upper signs and that with the lower signs is 
taken: 

= Vexpli [ dz^A + (x + ,x_ + z_)X p + j°{r J kj - m}p- x 
xVexplif dz + A^(x + + z + , x_ + \/2£o) i , 

Gf'Js = Pexpjz J Q dz + A^{x + + z + ,x^p- 1 °{ 1 J kj-m}p + x 
xVexplif dz_A + (x + + v^o,^- + -2-) 1 • 

[ JV2(yo-x -£o) J 



(34) 



The second order is sketched in figure EJ The kinks appear where the m^-terms are 
inserted. Now, there are not only contributions from on the backward light-cone but also 
from everywhere within. The different orders in the expansion are weighted by powers of 
tut according to the number of insertions and thus with the same weight wherever the kink 
is positioned. 

For large transverse mass my the adequate expansion would be given by an expression 
analogous to the weak-field case (section III B 1|) but without transverse components of the 
field. 



2. Gauge fields with dominant longitudinal and temporal gradients 

Taking the transverse components of the gauge field into account, too, but where the 
dependence on time and the longitudinal components is still more pronounced, the previous 
first-order expressions are modified to yield: 
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q(.i)+- _ -pexp/i / dz-A+(x+,x- + z_)X p + M T (x + , v^fo + xJ)p~ x 

xVexplif dz + A_(x + + z + , x^ + V^Co) i > 

_ pexp/i / ^ + A_(x + + 2 + ,x_)lp">lr(v / 2^o + ^+,^-)p + x 

xVexpli [ dz-A + (x + + V2£o,x~ + Z-)\ , (35) 

[ JV2(yo-xo-ta) J 

with A^r(3 ; + ) £-) = 1°[l J kj — m + , y J Aj(x + , £_)], the contribution from the transverse vec- 
tor components. All ingredients needed in order to put together the full expression in analogy 
to equation (J32j) can always be read off from the first-order equations, here and in the fol- 
lowing. The main difference with respect to the previous section is that now the trajectories 
that have a kink at a given point are no longer weighted uniformly with 7°{7 J /cj — m} but 
according to the value of ■y°{ , j J Aj(x + , x_) + ^ J kj — m} at the kink (x + ,x_) or, in higher 
order terms, the various kinks (see figure BJ). 

In case the transverse components of the field should be dominant compared to the 
temporal and longitudinal components one can start out with equation (|2*^|) and resum the 
terms involving the transverse gauge field: 

G(x ,y ,x 3 ,k 3 ,k T ) = 
= Vexp i d£ G Ao3=0 (x ,£ ,x 3 ,k 3 ,k T )U 3 (x - £oh°h 3 h + 7 • A(£ ,x 3 )] x 

- J Ha 

x U 3 (£ - x )G AQ/i =o{t,o,xo, x 3 , k 3 , k T ) G Ao3=0 (x ,y ,x 3 ,k 3 ,k T ), (36) 

with the definition for the four-vectors without transverse components = (v°, 0, 0, v 3 ) and 
where: 

(if 

G Ao 3 =o(xo, Vo, x 3 , k 3 , k T ) = Vexp{^= d£ x 

{ V2 JV2(y -x ) 

x (p + Mr(x+,x^ + £o)p~e + ^°d s (x) + p - M T {x+ + £ , X-)p + e~^° a3 ^) } .(37) 

The most important difference between the equations (J37j) and (f3T)|) is the fact, that 
a decomposition into two terms involving exponentials without Clifford matrices in the 
argument is not possible. However, because of the alternating occurrence of the projectors 
p ± , the expansion of this path-ordered exponential leads to two contributions per order: 
one with p + on the left-hand side and the other with a p~. See, for example, the second- 
order term: Together with the projectors, the translation operators appear in an alternating 
manner. This ensures that only contributions from on or within the backward light-cone 
are taken into account. That can be seen best after putting all translation operators to the 
right: 
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» o\ -» if® if® 

GAo3=o( x o,yo,x 3 ,k 3 ,k T ) = —7= I d£ ^= dr] x 

y I J\/2(y -x ) y 2 J V2(j/o-zo-£o) 

x p+A^r(x+, Co + x-) M T (x + + r]o, £ + x_)p-e-^"° %(a:) + 
+ p-M T ^o + x+, x_) -M T (£o + z+, + z_)p + e + ^ oa3 W] (38) 

The fact that there are no contributions to the propagator from outside the backward 
light-cone, consistent with relativity, also persists in the general case (see section 111 Al and 
figure EJ) . 

3. Gauge field of one time-like curvilinear coordinate 



The Green's functions obtained in [1| were equal to the retarded propagator only for 
a field depending on one time-like curvilinear coordinate n ■ x with n 2 > 0. Now one is 
in the position to give the retarded propagators for gauge fields as functions of one space- 
(n 2 < 0) or light-like (n 2 = 0) curvilinear coordinate, too. In order to investigate the different 
situations it is only necessary to look at one special choice for n in each case, because n 2 is 
a Lorentz scalar. For the time-like case n M = (1, 0, 0, 0) had been chosen, for the space-like 

= (0, 0, 0, 1), and for the light-like = (1, 0, 0, 

As already stated before, the time-like case is equivalent to the lowest-order expansion 
in the spatial gradient of equation (|25|). A reason for this particularly compact form can 
be read off figure |H The contributions to the propagator may be projected along straights 
of constant time, i.e., the same coordinate for which the retarded propagator's boundary 
conditions and hence the variable for the outer integrations are given. 



4- Gauge field of one space-like curvilinear coordinate 

If the gauge field is merely a function of a space-like curvilinear coordinate, say 
equations (|35|) can be rewritten as: 



VexpliV2 [ dz 3 A + (x 3 - z 3 )\ p + M T (x 3 - £o)p Vexp\iV2 [ dz 3 A„(x 3 + z 3 ) \ 
V J to > \ Jya-XQ J 



£(!)-+ = 

Vexp<i\/2 dz 3 A_(x 3 + z 3 )\ p~Mt(x 3 + £o)p + V exp <iV2 dz 3 A + (x 3 -z 3 ) 



(39) 



where the variable in the argument of the gauge field has been kept as A(x 3 ). As can be 
seen from this expression and also as indicated in figure El the field is integrated back and 
forth over the longitudinal coordinate while the phase factors are still functions of time [see 
equation (|33j)]. To the contrary, the respective Green's function derived for this case in 0] 
which is not a propagator, exhibits a different form, because a different boundary condition 
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had been imposed there. 

5. Gauge field of one light-like curvilinear coordinate 

For a gauge field being a function of the light-like coordinate x_, the equations 
simplify to: 



Here, one of the two path-ordered exponentials per addend becomes a simple phase. It is 
always the one, which is integrated along z + with z_ kept fixed. As the entire argument is 
only a function of z_ it can be taken out of the integral which no longer needs path-ordering 
and thus only gives a time difference. In figure El these phase factors are the contributions 
parallel to the direction of projection. Therefore, they are each mapped onto a point. The 
solution presented above is of the form relevant in a field theory with standard quantisation. 
The respective expression for the homogeneous solution of the Dirac equation in the case of 
retarded boundary conditions along the light cone was given in 0]. 

6. Longitudinally boost-invariant gauge fields 

A case of particular interest in ultrarelativistic heavy-ion collisions is that of a gauge field 
which is invariant under longitudinal Lorentz transformations, although this requirement 
often imposes further limitations on the configurations of the gauge field. Independent of 
this reasoning, usually also a dependence on the transverse coordinates is present. Here 
the simpler case is sufficient to clarify the main features. As in the forward light-cone 
r = \Jt 2 — z 1 > 0, it does not lead to ambiguities to replace it by half its square 
A± = x±a±(x + X-) and A T = A T (x + x_). Therefore the equations (J55|) become: 




It can be seen from these expressions that the longitudinal boost-invariance of the vector 




(40) 





(41) 
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potential does not lead to great simplifications. In fact the propagator itself does not be- 
come boost-invariant. The reason for this is that the boundary condition for the retarded 
propagator is not invariant under longitudinal boosts. From figure [7| it becomes evident as 
well that the problem does not exhibit the symmetry of the field. Nevertheless, there are 
simplifications for the homogeneous solution Gh- 



D. Ultrarelativistic heavy-ion collisions 



A standard way to describe an ultrarelativistic heavy-ion collision is to begin with a 
current of colour charges moving along the light-cone |l6[. In order to obtain the classical 
gauge field the Yang-Mills equations have to be solved in the presence of this current. In 
Lorenz gauge the field consists of two different contributions. 

One of them is made up of the Weizsacker- Williams sheets. These correspond to the 
Coulomb fields of the different colour charges boosted into the closure of the Lorentz group. 
Therewith they become 5-distributions in the light-cone coordinates. This approximation 
is justified for the description of particles with low longitudinal momentum (mid-rapidity). 
On the one hand they have a low resolution in this direction. On the other they cannot be 
comovers of the charges on the light-cone. For comoving particles the plane wave expansion 
would not be applicable 0. 

The other contribution is the radiation field in the forward light-cone A rad . By virtue of 
longitudinally boost-invariant boundary conditions on the light cone, it is sometimes taken 
as being invariant under longitudinal boosts, too. However, this is not necessarily the case. 
But anyway, as has been seen in section III C 61 there are no fundamental simplifications 
for the propagator. Thus this assumption will not be taken into account in the following 
discussion. 

In Lorenz gauge, the Weizsacker- Williams contribution to the gauge-field takes the form: 



A^ w (x) = -f E t a (t L a ) n j[x„-(b L _) nL 

n L =l 



A ww (x) 



g_ 

2tt 



E ta(ta)nj X + - (b 



In A 
In A 



n R =l 



(42) 



There is no contribution to the transverse components. A is an arbitrary constant, regular- 
ising the logarithm which does not appear in quantities like the field tensor. The t a are the 
generators of SU (3) c . The {t^ )R ) n represent the colour of the charges. 

In general, the presence of a second nucleus leads to a precession of the charges of the first 
nucleus and vice versa which is a manifestation of the covariant conservation of the current. 
This leads to further sheet-like contributions to the gauge field which can be combined with 
the Weizsacker- Williams fields by modifying the charges. Up to the next-to-leading order in 
perturbation theory they read |l6[: 
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(t L a)n L 



N R 

n,R=l 
N L 

(CZ + aS E fa b c{t L b )Zitc)nR d 
n L =l 



(bl) nR \n\\x T -(b«) nR \ + 0(a s 



-{b L ) nL \nX\x T -{b^) nL \ + 0{a s 2 ). 



(43) 



with the initial colours {t^ ,R )™ L . To all orders, the modifications amount to Wilson lines 
over the gauge field along the branches of the light cone. Hence, the precession terms are 
absent in adequate, i.e., light-cone gauges. In order to ensure colour neutrality of each 
nucleus the following relations have to be satisfied: Z) n £=i *a(*a)nl = = Sn^i^^a)^' 
Frequently, pairs of two charges belonging to the same nucleus are assigned to each other 
to form a colour-neutral dipole. Ultimately, higher-order perturbative calculations for the 
gauge field as solution of the Yang-Mills equations are only of value in the presence of 
hard energy- or short time-scales. Up to now, non-perturbative solutions for this problem 
have been obtained only in transverse lattice calculations |T3] • In any case, the general 
form of the field - continuous inside the forward light-cone and singular but integrable 
on the forward light-cone - allows to express the retarded propagator by a finite number 
of addends. Resumming the terms involving the Weizsacker- Williams-like contributions 
together with the longitudinal derivatives in the homogeneous part in equation (fTH|) yields: 



Pexp [i f ° dQf " ;- ><),{.,■) + ^kj + 7 • A ww (£ 0l x ) + 7 • A rad (£ , x) 
^ Jyo 

V exp i d£ V exp \i / dz j°[ij 3 d 3 (x) + 7 • A ww (z , x)} x 
x 7 [i 7 J 9j(x) + 7 j kj + 7 • A rad (£ 0l x)-m]x 
xVexpli / ? ° dz ^°[i^ 3 d 3 (x) + 7 • A ww (z ,x)} 



— m\ 



xVexp\i / d^\i 1 i dz{x)+ 1 -A ww {^x)\ 



(44) 



The path-ordered exponentials containing the Weizsacker- Williams contribution to the vec- 
torpotential can be written in the form of equation (jHOj) . because they do not contain any 
transverse components. Note that k% = has to be used as the longitudinal momentum 
has been kept inside the outer path-ordered exponential. Hence, after replacing A ww by 
equations (|4*2j) and carrying out the path-ordered integration, one obtains: 



( rxo 

Vexpli d^y^dsix) + 7 ■ A ww {^ x)} 

I Jyo 



'yo 

= p + V ] [ exp 

n L =l 
Nr 

+ p~v Yl ex p 

riR=l 



){o\ 



Xa 



(b L -)n L ] ~ 8[y- - (b L _) nL }}~ 
(&+)»«] -%+-(&+)»«]} 



-{xo-yo)dz(x) 



+ 



^+{xo-yo)d$(x) 



(45) 
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with 27ra^ = —gt a (t^' R ) nL R \n \ x T — (b T ,R ) nLR . It has been assumed that the sums in 

equation are ordered in such a way that, if m L > n L then (b^) mL > (b^) nL and likewise 
if m R > n R then {b R ) mR > (b R ) nR . The same is to hold for the products in the previous 
equation and has been indicated by the V. In contrast to equation (|3U|) . the translation 
operators are explicitely shown above. There they had been absorbed in equation (|28jl . 

Making repeated use of the relation f[9(+x)] = 9(+x)f(l) + 9(—x)f(0) and subsequent 
simplifications lead to: 



Pexp \i r d^[i^d,{x) + 7 ■ A ww {^ x)] 



' yo 

P + v n [i +{*[*-- 

n L =l 

x {%o - x ) exp[+ia^ L (x T )] 

N h 



)n B \ 

n R =l 



X 



%-- 








9(x - 


2/o)exp[ 


~i a n L 


(£r)] 


9[y+- 


" 






9{x Q - 


2/o) exp| 




(*t)] 



e -(sco-Vo)9s(a!) 1 



+ n _ [1 +{*[*+- $ 

e +(so-?/o)93(a;) /^g\ 



In ultrarelativistic heavy-ion collisions, the colliding nuclei are highly Lorentz-contracted. 
Therefore, in the infinite-energy limit 7 — > 00, the "longitudinal impact parameters" 
(b± R ) nRL go to zero like 7 _1 . However, care has to be taken, as important, because qualita- 
tively different contributions to observables could be neglected, although the expression for 
the propagator in this limit would still be correct. Taking the limit, the factors involving the 
differences of the step functions become the same for all terms of the first and the second 
contribution respectively, and thus can be factored out. Afterwards, the idempotency of the 
step functions remaining in the product can be used to split each contribution into three 
addends which can be recombined in pairs: 

Vexplt rd^ 1 °[t 1 %(x)+ 1 -A ww (^x)}}=Vexp\- f° dfryVW } + 

I Jyo ) I Jyo ) 

+ T L (x, ±)[0(x_) - %_)] e -(«°) a ^) + T R (x, ±)[6(x+) - 9(y + )]e + ^°~ yo ^ x \ 

(47) 



with: 



Nl 

T L (x, ±)+p + = P + (±1) Nl V J] exp[±m^(f T )] 



ul 

n L =l 



Nr 

T R (x, ±)+p~ = p-(±l) N *V I] exp[±i<(£ T )], (48) 



n R 1 

n R =l 



where ± in the argument indicates the sign of xq — yo- If the two nuclei are modeled by 
means of dipoles Nl and Nr are even and the overall factors with the sign are obsolete. 
What is left of the longitudinal structure of the nuclei is the path-ordering of the products. 
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It would become irrelevant after averaging over all permutations of the charges in which case 
it would be possible to write the product over the exponentials as simple exponential over 
the sum of the contributions. 

Putting equation (|T7)l into equation (pHj) allows for the reconstruction of the homogeneous 
part of the retarded propagators in the radiation field G r ^ d \ 

G™ d (x, y, k) = Pexp £° ^o7°[h j ^(x) + 7^ + 7 • A rad (£ , x) - m]] (49) 
and the final result reads: 

G H (x,y,k) = G^ d (x,y,k) + 

+ J d^G r H ad (x^,k)[T L ^,±)S(^) + T R (^±)5^ + )]G^,y,k) + 

+ J d^ v G r H ad (x, £, k)T L (i, ±)6(^)G r H ad (£ } 77, k)T R (rj, ±)5( V+ )G r H ad ( V} y, k) + 

+ J d^-qG^ix, i, k)T R {£, ±)5{i + )G r ^ d {i, 77, k)T L {r), ±)5(T]^)G 7 g d (r], y, k). (50) 

When the retarded (advanced) propagator is constructed from this homogeneous part, the 
minus (plus) sign in the argument of the T has to be taken. Higher orders in the T cannot 
contribute, because lines of constant x_ or x + can only be crossed once (see the discussion 
above). In the case of two colliding dipoles taking the ultrarelativistic limit in order to 
obtain equation (|4~7j) keeps only integration paths like the outer two depicted in figure |H1 
The other, qualitatively different paths, in which interactions with left- and right-movers 
are intertwined are neglected. With finite "longitudinal impact parameters", every charge 
would have to be treated on its own and the previous equation would contain contributions 
with up to Nl + Nji insertions. 

For a given field, the retarded propagator in the radiation field could now be obtained 
with the methods presented earlier in this section. The exact solution for the radiation field 
would be needed. Finally averaging over all configurations (current and hence field) would 
yield the desired observables. However, the transverse lattice calculation required at that 
point is beyond the scope of the present paper. 

III. FERMION-ANTIFERMION PAIR PRODUCTION 

In this section, the findings for the full propagator in a space-time dependent field are 
interpreted in view of the problem of particle production by vacuum polarisation. Such a 
calculation can be understood in two different ways. Either the field is an external field 
in the strict sense, or the field is the self-consistent solution of a set of equations. In the 
first case it is governed entirely by the characteristics of the physical system without the 
phenomenon of particle creation. Subsequently it can be used to determine the spectrum 
of produced particles. An approach of this kind is justified if particle production represents 
only a perturbation which must be checked eventually. In the second case, given the solution 
for the set of equations in form of the classical field is known, the question is how many 
particles were created in the process. Here, the set of equations would include the Yang- 
Mills equations with an initially present current and the current induced by the fermions and 
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antifermions. In this section the first way is to be pursued, because the back-reaction in the 
investigated system is suppressed by higher powers of the coupling constant not accompanied 
by powers of the vector potential. 

According to |l8| the expectation value for the number of produced pairs can be deter- 
mined by starting out with a negative energy fermion with four-momentum q corresponding 
to a positive energy antifermion: 



Mv) = v(q)e-^, (51) 

with the unit spinor v (q). After its propagation through the field, calculate its overlap with 
a positive energy fermion solution with four-momentum p: 



%{x) = u(p)e +i ^, (52) 

with the unit spinor u(p). Wave-function solutions of the Dirac equation at x and y are 
connected by the retarded propagator according to: 



tf;(x) = I d 3 yG R (x,yh%(y). (53) 
Hence, the required overlap reads: 



M(p, q) = J d 3 xd 3 y^ p {x)G R {x, y)-y%(y). (54) 
With the definition of the retarded one-particle scattering-operator T R 

G R (x, y) = G° R (x - y) + j d^d± V G R (x - t)T R (Z, V )G R ( V - y), (55) 
the overlap can be rewritten as: 

M{p,q) = J d^d 4 V M0^v)Mv) = u{p)T R (p,-q)v{q), (56) 

because ip q and ip p are orthogonal to each other and are solutions of the free Dirac equation, 
whence: 



d 3 yG R (x-yh My) = M*) 
J d 3 yMyh°G R (y -x) = $ 9 {x). (57) 

In the last step the definitions of ip q and ip p in equations ()51|) and (|52jl had been inserted and 
the integrations carried out. This led to the Fourier transformation of the single-particle 
scattering operator. The square of this overlap \M(p, q)\ 2 summed/averaged over all colours 
and spins is equal to the differential expectation value for the number of produced pairs 
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where the particle carries the momentum p and the antiparticle the momentum q. After 
integration over the phase space of the particle (antiparticle) one obtains the spectrum of 
produced antiparticles (particles). Integration over all momenta finally gives the expectation 
value for the number of produced pairs: 



where u p = \J\p\ 2 + m 2 and u q = y\q\ 2 + rn 2 . 

In order to address the issue of gauge invariance, rewrite the previous expression as 0: 

— <ftxd*ye? q '* c ~ v '\x {(7 • q + m)[ry ■ d(x) — m]G + -(x, y)[i^ ■ d(y) — m]} , 
(2~) lu.',, 

(59) 

where the derivatives are only acting on the propagator. The only coloured and hence gauge- 
dependent part of equation (J5T?j) is the propagator, even after taking the trace over its colour 
indices: 



tr {G(x, y)}^tr {G'(x, y)}=tr {Q(x)G(x, y)tf{y)} + tr {G(x, y)} (60) 
1. Gauge- dependent states / Wilson links 

Sometimes, the gauge invariance of the previous expressions is ensured by replacing the 
gauge-independent free states in equations (joTj) and by gauge dependent states ip, 
especially fermion and antifermion solutions in pure gauge fields, i.e., vacuum configurations 
of the gauge field: 

{i 7 • d(x) + Q(x) [7 • d(x)tf(x)] - m} i){x) = 0. (61) 

The gauge transformations associated with them cancel those coming from the propagator 
in equation exactly: 

M(p,q) = J d 3 xd 3 y^ p (x)G R (x,y) 7 %(y) - J d 3 xd 3 y^ p '(x)G' R (x, y)j%(y) = 

= J d 3 xd 3 y^ p (x)n\x)n(x)G R (x, y)rt(y)>fn(y)j> q (y) = M(p, ?) (62) 

A priori, the gauge in which if>(x) coincides with a free solution can be chosen arbitrarily, 
whence it must be specified in order to allow for an interpretation, which therefore is not 
unique. This goes hand-in-hand with the loss of an interpretation as states with a definite 
number of partons. In all but the one gauge of reference an arbitrary number of additional 
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gluons is included in them. Further, when deriving equation (|58|) q and p were the momenta 
of the fermion and the antifermion, respectively. Here, this correspondence remains obvious 
only in one gauge, which thus would be singled out. 

This method is equivalent to including Wilson lines with the correlators. For an ade- 
quate choice of these gauge links - especially their end-points - the gauge transformations 
connected to them cancel with those of the two-point functions. There the supplementary 
degree of freedom is the path along which it is evaluated. 



2. Average over the gauge group 

A way to guarantee a gauge invariant-expression without sacrificing the requirement 
of purely fermionic states is given by averaging the expression (|59|) for the expectation 
value over all gauges, because every supplementary gauge transformation u(x) can then be 
absorbed by shifting the functional-integral measure: 



[dft]tr {u{x)G{x, y)uj\y)} = J [d(wfl)]tr {G(x, y)} = I [dQ']ti {G{x, y)} . (63) 

As a first step, carry out a coordinate transformation from x and y to £ = x — y and 
C = (x + y)/2. Subsequently look at d(n) /<i 4 £. Any SU(N) gauge transformation can be 
parametrised according to 



N 



Q(x) = exp{-it a 6 a (x)} = H x )) (n(x)\exp{-i\ n (x)} 



(64) 



71=1 



with t a the generators of the gauge algebra and where X n (x) are real eigenvalues and \n(x)) 
orthonormalised eigenvectors of the matrix t a # a (x)^|. Functional integration of d (n) /o? 4 £ 
over the eigenvalues in equation (|60|) yields a non-zero result for x = y: 



d (n) 



d 4 £ 



d 3 q 



-d 4 ( x 



^ (2tt)3 2 ^ 

xtr{( 7 -g + m) [i 7 • <9(Ci)| Cl=c - rn] G?+_(Ci,C2) [* 7 • ^(Cs) l Ca =c ~ m ] } 



(65) 



Already in its present form the last expression is gauge invariant. Averaging over the eigen- 
vectors \n(x)) is not needed, because the trace already takes care of colour neutrality. Now 
the two gauge transformations in equation (JBOj) are evaluated at the same point and drop 
out of the trace. 

This result resembles calculations on the lattice. There the results contain the four- 
dimensional lattice-volume V as an overall factor. Hence, in that case one would look at a 
quantity like (n) /V. 

With the integrand in momentum space, the previous result becomes: 



d{n) 



d 4 £ 



d 3 q 



dH 



(27r) 3 2a; 9 (2vr) 



-tr 



{( 7 -q + m)G R -\l)G + 4l,l)G A -\l)} 



(66) 
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instead of 



(n) 



d 3 q 



(2tt) 3 2o;, 



-tr 



{(7 • q + m)G° R \q)G + 4q, q)G A \q)} , 



(67) 



which would be obtained directly from equation (jSHJ)- 

Making use of equation © (oi = +, a 2 = — ) as well as |H[2(J: 



r d 4 k 

T A (q,p) - T R (q,p) = J J7^ 5 ( k2 ~ m2 )^(k )T R (q, fc)(7 • k + m)T A (k,p) (68) 



and: 



gives: 



Gs{q,p) 



d 4 k 



5(fc 2 - m 2 )G^(g)T R (g, fe)( 7 • k + m)T A (k } p)G A (p), 



(69) 



<i (n) 



ci 4 g 
(2^ 



6(+q )2n5(q 2 - m 2 



d 4 k 



6(-k )2n5(k 2 - m 2 



d 4 l 



(2tt) 4 v uy v y (2vr) 4 
xtr{( 7 • g + m)T R (l, fc)( 7 ■ fc + m)T A {k, I)}. (70) 



This equation can be interpreted as the fermionic current induced by the external field from 
which the on-shell component (q 2 = m 2 ) is projected out with the extra condition that the 
current's momentum is on the mass shell once ik 2 = m 2 ). Opposed to that, in equation (J58|) 
the momentum is forced to be on-shell twice. 

Under a local gauge transformation Q(x), the full propagator transforms homogeneously 
[see equation (|60|)]. The considerations concerning gauge invariance remain unchanged, if 
this behaviour is preserved in an approximation. In the case where, in a given gauge and/or 
Lorentz frame a certain expansion should be exact, this requirement is trivially satisfied. In 
general, if the expansion is based on powers of the covariant derivatives and/or the mass, 
covariance is always ensured. This is the case for equation (|30|) . for expansions of equation 
(|32|) if the transverse momentum is omitted, and for approximations like equation (|35|) if 
the transverse derivatives are included with M. T [after adapting equations (|28|) and ([29)1] . 
If the approximate propagator should not transform homogeneously, still a gauge invariant 
result can be obtained by averaging over all gauges. However, in that case the situation has 
to be reinvestigated, because the condition x = y need not be sufficient to guarantee gauge 
invariance there. 



A. An example 

The radiation field in an ultrarelativistic heavy-ion collision has to be determined in 
transverse lattice calculations which are beyond the scope of the present paper. Hence, the 
previously determined expressions for the propagator and the quantities derived from it are 
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to be evaluated for a model for the radiation field. Along the lines of the discussion in 
section Hi Dl a longitudinally boost-invariant field is chosen: 



A ± = ± a ^e- x+x - /T ° 2 . (71) 
To 

The Lorenz d ■ A = and the Fock-Schwinger x ■ A = gauge-conditions are satisfied by this 
field. In its presence the one-particle scattering-operator defined in equation (f53J) with the 
propagator according to the lowest order expansion in the transverse components in section 
IHia reads: 



T R (p, q) = (2tt) 2 6®(p T - q T )T R (p, q) (72) 

with: 

V2T R (p, q) = 7 + T+(p, q) + l-T R (p, q) (73) 

where: 

rr± ( ~ ~x ^ f(r + r_T 2 ) 
t r [P, -Q) = Tar 



(r + r_r 2 )r± 

_ fl2 (jtorp)" (TiaroY /[r + r_r 2 /(^ + v + 2)] - f[ P± r T r 2 /Jjx + v + 2)] 
a V M! ^! r ± 2 [( yU + l)g T -(z/+l)p T ] 

(74) 

with r± = p± + g±. In equation (|72j) the omnipresent Dirac-distribution, originating from 
the transverse homogeneity of the field, is factored out. The expressions for the two compo- 
nents of the one-particle scattering-operator can be resummed partially into hypergeometric 
functions. / involves the exponential-integral function Ei: f(a) = a 2 e~ a Ei(a) — a. In ac- 
cordance with the discussion in the previous section, the approximation to the retarded 
propagator used here transforms like the full solutions, i.e., homogeneously. Higher order 
contributions to the retarded propagator G R (x,y) are also available analytically. They are 
qualitatively similar to the lowest order investigated here. However, the terms belonging to 
the different orders cannot be resummed easily. Further, a related approach to the problem 
of particle production investigated in [l| turned out to be in good agreement with the full 
solution. 



1. Gauge- dependent states / Wilson links 

An interpretation of equations (|72H74j) put into the integrand of equation (J58|) : 



\u 



(q)T R (q, -p)v(p)\ 2 = £ [8(p • T R )(p ■ T R f) - 4(p • q - m 2 ){T R ■ T R % (75) 
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along the lines of the discussion in section IIII II would be possible if the gauge of reference 
was chosen to be the Lorenz or the Fock-Schwinger gauge. Here, in these gauges, the vac- 
uum configuration of the gauge field vanishes identically per definition, whence the fermion 
solutions ip in the vacuum configuration coincide with the free solution. The last expression 



is summed over the eigenvalues A n of the colour matrix a = J2 n ^n\ n ) ( n \ 21 . 

In the case, where a light-cone gauge, say A + = was taken, one would have to resort to 
the expression of the one-particle scattering-operator in that gauge: 

Tr = 
and 

/(r + r_r 2 ) 



t r(p,-q) = -2ar - 



r + r_r 2 )r H 



_ 4fl2 v (+2mr )" (-2mT )" f[r + r_r 2 /(// + u + 2)] - /[ P+ r_r 2 /(/x + u + 2)] 
a y\ v\ r + 2 [(/i+l)g_-(z/ + l)p_] 

(76) 



In other words, the Wilson line taken along the ^--direction becomes unity. 



2. Average over the gauge group 

For a calculation according to section IIII 21 one can start in any gauge. In momen- 
tum space, the advanced single-particle scattering-operator is linked to the retarded via 
Ta(p, O) = ~ ^r(<2S£>)- Starting from equation (J75|) the integration over the momentum I and 
the trace can be executed to yield: 



d 2 l 



8q + k. 



tr{( 7 • q + m)T R (l, fc)( 7 • k + m)T A (k, I)} 
2iar ) u 



(2vr) 2 

+2iar Y 



x 



X 

fll v\ 

r] 2 (fi + 2) + 2ir 2 k_r] 
^2r 2 (fi + 2) 2 (fi + u + 3) 
j r\ 2 {y + 2) + 2iT 2 k-r] 
\2r 2 (z/ + 2) 2 (/x + i/ + 3) 



+ 



-4i(ar ) 3 | 

■ 4i(ar ) 3 | 
16(ar ) 4 

/i + z/ + 4 1 4r 2 (z/ + 3)(/i + 3) 



-r 2 A;_ 2 m 



-r 2 kJ\n 



l + i 



^ + 2) 
k-r 2 

yjis + 2) 
kW 



r\ 2 [y + /i + 6) 4ik^rj(p + fi + — fi) 



+ 



2(z/ + 3) 2 (/x + 3) 2 



r 2 k. 



3) : 



In 



.7^ + 3)" 

' A;_T 2 



TQ 2 k_ 2 

(/i + 3) : 



In 



l + i 



k-r 2 



+ 



(77) 



The quantity r\ regularises the expression. The divergencies are due to the symmetry (boost 
invariance) of the field and to the expansion of the propagator in the transverse mass. It 
can be checked - e.g. by interchanging the indices v and \i - that the last expression is 
in fact real, as it should be. Further the sums correspond to hypergeometric functions and 



24 



their special cases like gamma functions. The logarithms of complex argument are linked 
to logarithms of the modulus of the argument or arcus tangents respectively. However the 
above is still the most compact representation. The final result has to be summed over the 
eigenvalues of the colour matrix a as before. 

The momentum k cannot be identified with the antiparticle momentum p. This can be 
seen after integrating over the plus and minus components of the momenta k and q and 



taking the resulting expression at vanishing transverse mass y \kx\ 2 + m? = 0: 



(2tt)V d (n) 



L t 2 t] 2 k 2 d^d 2 k T d 2 q T 



(arof + E (+2mT ° r (_2mT( 



£,=0=^J\k T \ 2 +m? V V- V - 

( 2t{aT ) 3 {u - fj,) 4(ar )> + /i + 6) 



(ji + 2)(v + 2)(ji + u + 3) (y + 3)(ji + 3)(ji + v + A)) ' 

(7s; 



Use has been made of the relation (2ir) 2 [8^ 2 \kT — It)} 2 = LT 2 S^ 2 \kT — It) with the trans- 
verse length Lt- k serves to regularise the integrals. In the final result r\ and k occur 
exclusively in the dimensionless product rjK. Due to Noether's theorem, in the absence of 
transverse components and variations of the field, the particle and antiparticle should be 
produced in strict back-to-back configurations only. Hence, the right-hand side of the last 
equation should have support merely for q T = —k T = 0, but this is not the case. In fact it 
does not depend on qT at all. Already in the general expression (ffO|) . the dependence on the 
particle momentum q is only via the phase-space volume and the factor in the trace. There- 
fore, this quantity can only be interpreted consistently after integration over the momentum 
variable k. For the present example, all the integrals could be carried out analytical. Here 
the result for massless fermions in a system with one temporal and one spatial (longitudinal) 
dimension is spelled out, where it is the exact result: 



(27r)V d (n) 



7] 2 K 2 



d 2 £ 



(1+1) 



4(ar ) 2 + ^[l 



cos(2ar )] 



(79) 



For strong fields and/or slow decay-times, the gauge-field (fTTj) corresponds to an (almost) 
constant electrical field. In that limit the last expression is dominated by the quadratic 
term. This observation is in correspondance with the behaviour of the constant-field result 
by Schwinger 0] in the case of vanishing transverse mass. 



IV. SUMMARY 

The creation of fermion-antifermion pairs due to vacuum polarisation in the presence 
of classical fields has been investigated starting from the advanced and retarded fermion- 
propagators in the field. The needed propagators were derived as Green's function solutions 
of the Dirac equation to all orders in the arbitrarily space-time dependent field. The re- 
quired propagators were chosen by imposing the respective boundary conditions. In the 
course of the calculations, [/(A)-generators have been admitted for the charge-space of the 
gauge field. Thus QED and QCD effects can be taken into account simultaneously. In terms 
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of creation processes electron-positron production as well as quark-antiquark chromo- and 
photo-production can be adressed. After a discussion of the general features of the propa- 
gator and its constituents, like the generalised translation operator, its characteristics have 
been treated in great detail. In order to gain further insight into the structure of the ob- 
tained Green's functions and in order to get the link to the results for purely time- dependent 
fields in jlj , various expansion schemes have been derived from the full expression. Namely, 
these are the weak-field, the strong field, the gradient, and the Abelian approach as well 
as some combinations thereof. Further, also for the same reasons as above, the features of 
the propagator are investigated thoroughly in different classes of field configurations. This 
includes situations that are dominated by temporal and longitudinal degrees of freedom and 
last but not least longitudinal boost-invariance of the system. This feature is frequently of 
importance for the description of ultrarelativistic heavy-ion collisions in the semiclassical 
limit. Therefore this case has been spread out in detail in a seperate subsection. How- 
ever, it has been found that the assumption of longitudinal boost-invariance of the gauge 
field does not lead to significant simplifications for the propagator. The reason herefore is 
that its boundary condition does not display this symmetry. Contrary to that there are 
simplifications for the homogeneous solutions. 

The issue of the gauge invariance of the expression for the expectation value for the 
number of fermion-antifermion pairs produced in the presence of the field is discussed in 
detail. These expressions are functionals of the propagators which have been derived above. 
These two-point functions are the gauge dependent (coloured) objects. Even the trace over 
them is in general not gauge invariant. Hence, two different ways to generate gauge-invariant 
quantities are discussed. On the one hand, this can be achieved by projecting onto gauge- 
dependent states instead of free fermion solutions. This method is equivalent to the inclusion 
of Wilson lines. In a given gauge of reference the new expression coincides with the bare 
one. This is where it is evaluated. Nevertheless, the interpretation as purely fermionic 
states is lost, because in all but the gauge of reference, in the assymptotic states gluons 
are included with the fermions. On the other hand, in order to preserve purely fermionic 
states, one can average the expectation value over the gauge group. In both approaches, it 
seems not to be possible to identify the particle- and antiparticle-momentum simultaneously. 
While the momentum of the particle can be identified, in the second case, the supplementary 
momenta, do not satisfy Noether's theorem. In the first case, the additional gluons turn 
the momentum-like variable more into a continuous parameter for the numbering the states 
than into the antiparticle momentum. 

Finally, the various expressions are evaluated in the presence of a model for the radiation 
field in an ultrarelativistic heavy-ion collision, because the exact field configuration, for 
example in the McLerran-Venugopalan model is only availabe numerically. For said field, an 
analytic expression for the one-particle scattering operator to lowest order in an expansion in 
the transverse components of the covariant derivatives and the mass is given. This expansion 
preserves the behaviour under gauge transformations of the considered quantities. Higher 
order terms for the propagator in this field are also available analytically. In the framework 
involving the gauge-dependent states (Wilson lines), after carrying out the traces over the 
colour and the Clifford algebra, this result provides a two-parameter inclusive distribution 
of the states as explained before. For the method, where the average over the gauge group 
has been taken, the integration over all momentum variables can be carried out analytically. 
In a certain limit, the present model field-tensor becomes constant. There, a behaviour 
quadratic in the field strenght is obtained. For vanishing transverse mass, the same holds 
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for the constant-field Schwinger-f or inula. 
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FIG. 1: Contour of integration in the complex time-plane for the Schwinger-Keldysch formalism 
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FIG. 2: The two contributions to equation (IM()(1 from the backward light-cone. The white line 
corresponds to the path-ordered exponential along the black to that along £ + . The origin of 
the coordinate system is situated at (xq,xs). The dashed line represents t = yo- 
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FIG. 3: The two contributions to the equations or (|35|) respectively but for the second order. 
The white lines correspond to path-ordered exponentials along Z-, the black to those along z + . 
This figure also represents equation IjUHJ) but where the black and the white lines stand for simple 
integration-paths. In all cases, the kinks are found at t = £o an d t = 770, i.e., where the "transverse" 
matrices are inserted. The origin of the coordinate system is situated at (xq,xs). The dashed lines 
represent lines of constant time. The dotted lines stand for the backward light-cone. 
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FIG. 4: The two contributions in the case of a purely time-dependent field. The sloped white 
lines correspond to path-ordered exponentials along z~, the black to those along z + . The kinks 
are found at t = £o and t = rjQ, i.e., where the Mt are inserted. The origin of the coordinate 
system is situated at {xq,x^). The dashed lines represent lines of constant time. Along these the 
black and white lines in the direction of the light-cone coordinates are projected onto the vertical 
ones at the sides of the figure. These projections represent the paths along which the path-ordered 
exponentials are actually evaluated in the purely time-dependent case. The dotted lines stand for 
the backward light-cone. 
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FIG. 5: The two contributions in equations (f.'fflj) but for the second order. The sloped white lines 
correspond to path-ordered exponentials along z~, the black to those along z + . The kinks are 
found at t = and t = rjo, i.e., where the A4t are inserted. The origin of the coordinate system is 
situated at (xo,^)- The dashed lines represent lines of constant time. The dotted lines stand for 
the backward light-cone. The dash-dotted lines mark straights of equal longitudinal coordinate. 
Along those the black and white lines in the direction of the light-cone coordinates are projected 
onto the horizontal ones at the bottom of the figure. These represent the paths along which the 
path-ordered exponentials in the equations (|39|) are actually evaluated. 
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FIG. 6: The two contributions in the equations (|40f) but for the second order. Of the lines on and 
inside the backward light-cone (dotted) the white correspond to path-ordered exponentials along 
Z-, the black to those along z + . The kinks are found at t = £o and t = rjo, i.e., where the Mt are 
inserted. The origin of the coordinate system is situated at (2:0,23). The dashed lines represent 
lines of constant time. The dash-dotted lines mark straights of constant z_. Along those the 
black and white lines inside and on the backward light-cone are projected. These projected lines 
represent the paths along which the path-ordered exponentials in the equations (|4Uj) are actually 
evaluated. Especially the black lines are projected onto points which leads to simple phases as the 
field remains constant. 
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FIG. 7: The two contributions in the equations (|41|) but for the second order. The white lines 
correspond to path-ordered exponentials along z—, the black to those along z + . The kinks are 
found at t = £0 and t = 770, i.e., where the Mt are inserted. The origin of the coordinate system 
is situated at {xq,x^). The dashed lines represent lines of constant time. The dotted lines stand 
for the backward light-cone. The dash-dotted hyperbolas are characterised by constant proper 
time r and hence constant values of the gauge field A. Mark that in general there can be an offset 
between the origin of the coordinate system for these curves (0,0) and that in which the propagator 
is obtained {xq,x^). Therefore, the propagator does not exhibit the symmetry properties of the 
field. 
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FIG. 8: The six contributions to the retarded propagator of highest order in the Weizsacker- 
Williams fields of two colliding dipoles. The solid lines depict the trajectories of the charges, 
the dashed generic integration paths which, in the forward light-cone, correspond to retarded 
propagators in the radiation field. 
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